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INTRODUCTION 


It has been shown in recent papers by Vohr and Chow (^Ref . if], Hirs f^Ref . 2^ 
and Malanoski 0Ref \ 3^ that spiral-grooved, journal bearings have perform- 
ance characteristics which are significantly better than those possessed by 
plain journal bearings. A particularly important asset of gas- lubricated 
spiral-grooved bearings is that they are capable of high stiffness and high 
load capacity at high A operation are due to the "self pressur izing"ef fee t 
of the spiral grooves. 


Because the differential equation derived for the pressure distribution around 
a spiral- grooved bearing is such a complicated one, the only solutions to the 
equation obtained to date have been based on perturbation analysis, and are 
valid for small radial displacements only. These solutions have proved valu- 
able in providing basic design information on spiral-grooved bearings, and 
because of their relative computation ease, will probably continue to be exten 
sively used to obtain optimization design data. However, in using small dis- 
placements solutions, the question necessarily arises as to how accurately 
they apply to bearings operating at moderate to high eccentricity ratios, A 
second question that arises is whether the values of groove parameters which 
produce optimum performance at low eccentricities will also provide optimum 
performance at large displacements. 

In order to resolve these and other questions, a numerical solution for the 
pressure distribution in a spiral- grooved bearing was obtained valid for 
arbitrary displacement and misalignment of the journal. Results of that solu- 
tion are presented in this paper. These include performance charts, giving 
the load capacity and attitude angle characteristics for one optimized con- 
figuration of grooved bearing as a function of eccentricity ratio and bearing 
number A. 

Other aspects of grooved bearings on which data are given are: the effect of 

length to diameter ratio on bearing load and attitude angle, the restoring 
moment generated by misalignment of the journal, the variation of optimum 






values of groove parameters with eccentricity ratio and the effect of groove 
depth on the linearity of the load vs. eccentricity curve. 

Present theoretical results are compared with recent experimental data of 
Malanoski . 



THEORETICAL AND NUMERICAL TREATMENT 


The data presented in this paper are the result of a numerical solution of an 
approximate form of Reynolds equation applicable to spiral groove geometries 
with a large number of grooves. This equation was established by J . H. Vohr 
and C. Y. Chow QRef. l] and has as the dependent variable the "smoothed 1 ' pres- 
sure distribution. The basic approximation made in the derivation is that the 
actual pressure distribution deviates from the smoothed one by some small amount 
over each ridge-groove pair and that the detailed shape of the deviation becomes 
unimportant as the number of grooves approaches infinity,, The detailed pressure 
gradient components over each of the ridges and grooves are considered constant 
so that the flow can be calculated and related to the smoothed pressure profile. 

With the nomenclature introduced by Fig. 1, the smoothed pressure equation as- 
sumes the following form: 
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where the M's are the flow rates in the 9 and z directions in the ridges and 
grooves; they are related to the pressure distribution by the following rela- 
tions : 
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In the present treatment, the speeds U and V and tnQ clearances h and h were 

r § 

specialized to the following forms: 
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where A - groove depth = constant. 

The values of x and y represent the coordinates of the shaft center at z = 0 
o o 

measured with respect to the concentric position. The angles a^, and a 2 allow 
misalignment of the shaft and bearing axes. The geometry is presented in Fig. 2. 

The numerical treatment of Eq. (1) follows the method introduced in references 
5 and 6. 


Eq. (1) is a second order partial differential equation with variable coeffici- 
ents and contains the first derivative of density with respect to time. The 
difference between this equation and all other forms of Reynolds equation en- 
countered in the literature is in the presence of the cross, second space 


derivative 


dzc>e 


The assumption that the film is isothermal leads to the replacement of the den- 
sity by the pressure P throughout. The equation is now non-linear, containing 
2 

P as well as P. 


The change in variables 

Q = P 2 (10) 

leads to the same form of equation in Q with the quantity V~Q in some of the 
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I 

coefficients. ‘ 

I 

Approximating all space derivatives of Q by standard central difference formulae 
and the time derivative by a forward difference formula, equation (1) becomes | 

algebraic, involving, at each point of a finite cell mesh drawn on the gas film, j 

the value of Q at the point and its eight immediate neighbors . 

i 

Letting the value of at each grid point be a part of the coefficients, j 

(imagining that it is known) the resulting algebraic equations are linear in 

Q. 


Due to the fact that none of the central differences introduce points further 
away from the center point than one column and one row, equation (1) reduces 
to the form 


Wj (Q}j + Hj + Hj (Q > j+1 = ( R )j ! j = '= N 


(11) 


where, for each value of j, the M equations pertaining to an entire column of 
points (i = 1, 2, . , .M) are written at once. Therefore, each matrix equation 
(11) contains M equations thus making the dimensions of > (jGj J GGj 

M x M and the dimensions of ^ 


be M. 


It should be noted that equation (11) can also be used for boundary points as 
long as the proper values axe inserted in the coefficient matrices and the 
right hand side. This is shown below. 


i 

I 


The boundary conditions are of the following types: 

a) given pressure (such as ambient) at the bearing ends. 



( 12 ) 


Equation (12) is of the same general type as equation (11) so that this 
type of boundary condition can be formulated together with the field 
equations . 
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b) flow matching between two grooved regions of different characteristics, 
This includes the case of grooved regions adjacent to plain regions due to 
the fact that, if equation (1) reduces to the equation valid for 

plain bearings. 


The present treatment is limited to cases where dissimilar grooved regions 

are joined on lines of constant z. Then the conditions can be expressed as 
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where + and - mean the two sides of the joint line The derivatives of 
pressure implied by (13) on either side of the joint line are written by 
means of two point formulae so that also equation (13) involves the joint 
point and its immediate neighbors and can be written in the same form as 
equation (11) . 


c) lines of symmetry. These only occur at z - constant (center plane of 
bearing) . This condition means that no flow in the z direction can exist 
at any 9 position or 

(M ) a + (M ) (1 - a) = 0 04) 

z g z r 

which can be used by writing the regular field equation at the symmetry 
point and eliminating all values of Q at the points on the non- treated 
side of the boundary by means of equation (14) 


d) cyclic condition 
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Letting j vary in the 9 direction, equations (15) and (16) reduce to the 
following: 
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However, equations (17) and (18) introduce two difficulties: first, they 

do not relate three or less adjacent columns of unknowns, and second, they 
introduce (q}^ and ($)q which are outside of the grid of unknowns. Refer- 
ence [0 suggests how to write the regular field equations on columns 1 
and N and gives a procedure for satisfying the boundary conditions (17) and 
(18). 

Therefore, the columnwise influence coefficient methods introduced in Reference 
can row be used to solve equation (11) . 

Two solution strategies can be followed: 

a) Iteration method. The time derivative is set to zero, the value of 

~\J Q in the coefficients is started at some guess, the value of Q is found 
at each grid point by solving equation (11), this value is used in V Q 
and the answer is obtained by iteration. This method is quite efficient 
when it converges because it usually gains one digit accuracy every one or 
two iterations. However, for severe conditions of eccentricity and bear- 
ing parameter A, severe numerical instability hinders the success of this 
method of solution. 

b) Diffusion net. hod: the time derivative is left in equation (1), an 

/’ 

initial distribution of Q is selected and used to evaluate the "V Q terms, 
and equation (11) is then solved repeatedly to find Q (AT), Q (2AT) , . , . 
where AT represents a suitable time step. Therefore, the progress of the 
pressure distribution in time is followed as in the natural diffusion 
process by which the pressure adjusts to a change in running conditions . The 
asymptotic value of the pressure at any point is that corresponding to 
steady- state , This method can always be made numerically stable by choice 
of a suitably small value of AT. Therefore, although more computing time 
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is necessary, answers can be gotten even for severe running conditions. 

The output data include the integrals of the pressure distribution so that 
the resultant forces along the x and y axes and the resultant moments about 


axes 

are known. 
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The computer program used for implement, a tion of the solution scheme is written 
in FORTRAN IV, and can accomodate slider bearings as well as journal bearings. 


RESULTS 


Steady-state , radial displacement , performance data calculated for the spiral- 
grooved journal bearing are shown in Figs. 3 and 4. Fig. 3 presents the 
dimensionless load W = vs . bearing number A = (^) for various 

values of eccentricity ratio e while Fig. 4 presents data on the bearing atti- 
tude angle ^ as a function of A and e. The data presented are for a bearing 
of L/D = 1.0, with grooved member rotating, having the following values for 
the groove parameters. 


p = 25° 

a ~ 0.35 (optimum groove parameters for maximum radial 

r = 2.33 load at A = 20) 

Z - 0.5 


These are the optimum values to provide maximum radial stiffness* at A = 20 
and € = 0 for the case of grooved member rotating. This particular bearing 
was chosen for high eccentricity analysis because it had been studied experi- 
mentally [Ref. 3.] and a stability chart had been prepared for it. It should 
be noted that optimum groove geometry for maximum stiffness at other values 
of A are not too different from those presented above. Optimum values for 
groove parameters are somewhat different if the smooth member rotates rather 
than the grooved member. In this paper, all data reported are for the case 
of grooved member rotating. 

In Fig. 5 is compared the load vs. eccentricity and attitude angle character- 
istics of a Pi ain bearing and our spiral- grooved bearing at A= 9.0. The 
plain bearing data were obtained from a study by Raimondi Q Ref . ?3 As can be 
seen, the non-linear increase of load with eccentricity for the optimized 
grooved bearing is very nearly the same as for the plain bearing. Also, the 


* Radial stiffness at e = 0 is defined as 


W cos 6 


where e < < 1 . 
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attitude angles for both bearings decrease in the same manner with e- One can 
note, however, that the attitude angle for the grooved bearing is lower than that 
for the plain bearing. 

When run at other values of A(i.e., A = 0.1 and A= 1.5) the load capacity of 
the grooved bearing was found to increase more rapidly with eccentricity ratio 
than does that of a plain bearing. As we shall see later, the linearity of the 
load vs. eccentricity ratio curve for spiral grooved bearings does depend on the 
values of the groove parameters. Hence, it is difficult to generalize about the 
linearity of the load curve other than to say that like plain bearings, the ra- 
dial stiffness of grooved bearings does always appear to increase with eccentri- 
city ratio. 

In Fig. 6 is shown the effect of length to diameter ratio on the performance of a 
grooved bearing at A = 1.5. Increasing L/D from 1.0 to 2.0 for the grooved bearing, 
results in an increase in dimensionless load of 90%. In comparison, the increase 
for a plain bearing would be only 75%. The greater percentage increase obtained 
for the grooved bearing is probably due to its self-pressurizing effect, which 
would be greater as length increases. One would expect this effect to be more 
significant at higher values of A and to be insignificant as A— >0. 

In the present study, the moment M that would develop on the journal due to mis- 
alignment was calculated for one value of A(A ” 1.5). The results are shown in 
Fig. 7. The abscissa in this figure is 5L/C where 6 is the angle of misalignment 
in radians. 6L/C represents the change in eccentricity ratio from one end of the 
bearing to the other. Also shown in Fig. 7 are misalignment moments for a plain 
bearing obtained from a linear perturbation analysis by Ausman Ref. 8 

There are two interesting features about the data in Fig. 1, One is that the 
misalignment moments for the grooved bearing are very linear with b (the moments 
plotted for the plain bearing are necessarily linear because they were calculated 
from perturbation analysis). The second is that the attitude angle of the mis- 
alignment moment is very much less for the grooved bearing than for the plain 
bearing (39° for the grooved bearing vs. 80° for the plain bearing). Also, the 
grooved bearing seems to have a larger misalignment moment relative to its load 
capacity than does the plain bearing. For example, th% ratio of the moment, M/L 
at 6L/C =0.1 to the load W at g = 0.1 is .0263 for the grooved bearing but only 


.022 for the plain bearing. These ratios can be used to estimate misalignment 
moments from load data for bearings with L/D *■ 1.0 if misalignment data is 
unavailable . 

An important question concerning the performance of grooved bearings at high 
eccentricity is whether the values of groove parameters that yield maximum 
radial load at low eccentricities will also give maximum radial load at high 
eccentricity. To investigate this, we determined the groove parameters to 
give maximum radial load for e < < 1 at A= 0.1, and then studied the effect 
that changes in these parameters would make on radial load at e = 0,6. The 
results are presented in Fig. 8a through 8d . The optimum values for (3, a, F 
and Z for e — £> 0 are indicated, respectively, by the dashed vertical lines 
in Figs. 8a through 8d . These values are: 


0 = 29° 

a -■ 0.51 (optimum groove parameters for maximum radial 

r « 2.2 load at A « 0.1) 

Z = 0.765 


As can be seen in Fig. 8, the above optimum values of 0, a and Z also are 
very close to the optimum values for e = 0.6. An exception is the parameter 
r = (h^/h^) £ _ q, the optimum value for which is approximately 1.8 for £ = 0.6 
That the optimum value of T should decrease with eccentricity is quite reason- 


able physically. It follows from the fact that to maximize bearing load at 
high e, one would want h^/h^ to at optimum value in the narrow clearance 
region where the load carrying pressure profile is developed. To keep h /h 


g‘ r 

at a local optimum value in a region where h and h are becoming smaller, 

S ^ 

one would have to reduce h - h i.e. reduce F = (h /h ) n 

g r g r' e « 0. 


If one used a value of 1.8 for JT in the design of a grooved bearing at A= 0*1, 
the radial load at low eccentricity would be reduced below that obtained with 
F = 2.2 whereas the load at e = 0.6 would be increased. Hence, the load vs. 
eccentricity curve would be even less linear than it would be with F - 2.2. 

The reverse possibility also exists. That is, by using a value of T greater 
» than 2.2 one should reduce load capacity more at high e than at low e and 


produce a more linear load vs. e curve. This performance characteristic would 
be desirable in many bearing applications e.g. gyroscope suspension systems. 


To try to improve the linearity of the load vs. e curve, calculations were 
made using r * 2.8 and 3.2. As an index of linearity we can form the ratio 
(M > 6 = 0.6 


(W >e = 0.3 


With perfect linearity, this ratio would be 2. From our calcu- 


lations we obtain that this ratio is 2.3 for F = 2.8 and 2.18 for F= 3.2 < 

For r - 2.2, the ratio is 2.57. We thus find that increasing F does signifi- 
cantly improve the linearity of the load vs. e curve. 


The above example points up one of the important advantages of spiral-grooved 
bearings, namely, their flexibility to be designed to satisfy many different 
performance requirements by proper choice of the groove parameters. Of course, 
this is not without its drawbacks. Due to the many design variables associ- 
ated with spiral-grooved bearings, proper design of these bearings can be a 
very complicated task. 

Before leaving the discussion of optimum values for (3, a, T and Z, we can note 
that the optimum values of these parameters for A= 0.1 are not too different 
from the optimum values for A= 20. More significantly, the bearing optimized 
for radial load at A = 20, when run at A = 0.1, has only 137» less radial load 
capacity than the bearing optimized specifically for operation at A = 0,1 
Hence this A = 20 design can give satisfactory performance over a wide range 
of A. 

As noted earlier, experimental data have been obtained by Malanoski [_Ref. ^3 
for the spiral- grooved bearing optimized for maximum stiffness at A = 20. 

The load vs. eccentricity data are shown in Fig. 9 along with the appropriate 
theoretical curves obtained from Fig. 3. Two test bearings were used, each 
of length and diameter equal to 1.5 inches with 36 spiral grooves etched onto 
the journal surface. 


As can be seen, agreement between experiment and the < is quite good at all 


r ■ 


r 


values of A up to 15.5, with the unexplained exception of the data taken at 
A= 2.68, which is quite a bit lower than theory predicts. The fact that all 
data lie somewhat below the theoretical curves can be partly explained by the 
so called ,? edge effect 11 considered by Muid jerman QRef . for incompressible 
thrust bearings. The predicted degradation of load due to edge effect in 
Malanoski' s bearing, as calculated from [Ref. lcT), would be 87„. 

In Table 1 below the values of attitude angle measured by Malanoski at higher 
values of eccentricity are compared with theoretical values obtained from the 
curves . in Fig. 4 . 


TABLE 1 

Comparison of Predicted and Measured 
Attitude Angles for Spiral Grooved Bearing 
L/D = 1.0, p = 25°, a = 0.35, T = 2.33, Z = 0,5 


Attitude Angle 


A 

e 

Expt 

Theory 

1.25 

.396 

58° 

52° 

1.25 

.523 

58° 

47° 

2.68 

.327 

50.5° 

41.5° 

2.68 

.487 

48.5° 

36° 

8.47 

.269 

21.6° 

17° 

8.47 

.310 

23.5° 

16.5° 

15.5 

.27 

11.2° 

7° 

15.5 

.333 

9.1° 

6.5° 


In general, agreement between theory and measurement is reasonably good 
although measured values of attitude angle are consistently larger than 
predicted . 


For sake of completeness, we have included in this paper a stability map 
for the spiral-grooved bearing tested by Malanoski. This map is 
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shown in Fig. 10. Discussion of the theoretical basis for this map is beyond 
the scope of this present paper; the interested reader is referred to the 
papers by Pan^Ref. 4 J and by Malanoski JjRef. 3^j for details of the analysis 
leading to construction of the map. 

The ordinates for Fig. 10 are bearing number A and bearing dimensionless criti- 
cal mass M . Use of the map for design purposes is very simple: one simply 

c 

determines the bearing number A at which one plans to operate, and then calcu- 
lates the maximum value of critical mass corresponding to that A. Actual rotor 
mass must be less than this value for stability. 

Alternatively, one may enter Fig. 10 with a dimensionless critical mass, and 
determine the maximum value of A at which the bearing can operate stably. 

The stability data in Fig. 10 pertains to unloaded bearings operating at e - 0, 
Limited experimental data indicates that loading the bearing improves its 
stability. Hlef. 111. 

Full experimental verification of the stability map in Fig. 10 has not been 
achieved. However, the test bearing of Malanoski was run unloaded at 60,000 
RPM(A = 20) with no sign of instability. The rotor mass was 2.125 lbs. per 
bearing which gave a value for the dimensionless mass parameter of 0.2. Pre- 
dicted critical speed for the bearing was approximately 80,000 RPM. 


SUMMARY AND CONCLUSIONS 


A solution for the pressure distribution around a spiral-grooved journal bear 
ing operating at large eccentricity ratios is obtained numerically by the 
method of influence coefficients. Performance charts giving load capacity 
and attitude angle are provided for a bearing optimized for radial stiffness 
at A 55 20 for the case of grooved member rotating. The charts cover the 
range of A from 0.1 to 20. Curves are given for eccentricity ratios up to 
e =0.7. The bearing considered has good load carrying capabilities over the 
entire range covered by the performance charts and is therefore a useful de- 
sign. A stability map is provided for the bearing. 

It is found that the load capacity of a spiral-grooved journal bearing optim- 
ized for maximum radial stiffness increases non-linearly with eccentricity 
ratio in qualitatively the same fashion as does the load capacity of a plain 
bearing. Also, the attitude angle for grooved bearings decreases with ac- 
centricity in essentially the same manner as for plain bearings. Hence, it 
is tentatively recommended that in cases where only low eccentricity data is 
available on a grooved bearing, this data can be extrapolated to high eccen- 
tricity by assuming that the grooved bearing behaves similarly to a plain 
bearing. Calculations indicate that this would result in a conservative es- 
timate of load capacity for the grooved bearing. 

The variation of grooved bearing performance with length to diameter ratio is 
briefly examined. If appears that the increase in dimensionless load with 
L/D is greater for grooved bearings than for plain bearings. 

The restoring moment generated due to misalignment of a grooved journal is 
found to increase approxima tely linearly with the angle of misalignment. The 
attitude angle of the restoring moment is much lower for grooved bearings tha 
the plain bearing. 

Optimum values of the grooved parameters (3, Q:, Z, to produce maximum radial 
stiffness at high eccentricity were found to be approximately the same as 
those obtained from small eccentricity analysis. However, it was found that 
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the optimum groove depth ratio, T, for maximum load at high eccentricity should 
be smaller than that producing maximum load at low eccentricity. It was also 
found that the choice of T significantly affected the linearity of the load vs. 
eccentricity, with larger values of T tending to make the curve more linear 
although with some sacrifice in load capacity. In some applications, however, 
linearity of bearing response may be more important than large load capacity. 
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